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ABSTRACT 

The cosmic star formation rate density (CSFRD) has been observationally investigated out to redshift z — 10. 
However, most of theoretical models for galaxy formation underpredict the CSFRD at z > 1. Since the theo- 
retical models reproduce the observed luminosity function (LF), luminosity densities (LDs), and stellar mass 
density at each redshift, this inconsistency is not simply implied that theoretical models should incorporate 
some missing unknown physical processes in galaxy formation. Therefore, we examine the cause of this in- 
consistency in UV wavelengths by using a mock catalog of galaxies generated by a semi-analytic model of 
galaxy formation. We find that this inconsistency is due to two observational uncertainties: dust obscuration 
correction and conversion from UV luminosity to star formation rate (SFR). The methods of obscuration correc- 
tion and SFR conversion used in observational studies result in the overestimation of CSFRD by ~ 0. 1 -0.3 dex 
and ~ 0.1-0.2 dex, respectively, compared to the one directly obtained by our mock catalog. We present new 
empirical calibrations for dust extinction and conversion from observed UV LF and LD into CSFRD. 
Subject headings: galaxies: evolution — galaxies: formation — methods: numerical 
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1. INTRODUCTION 

The cosmic star formation rate density (CSFRD), p*, is one 
of the most fundamental quantities that reveals how present 
galaxies are formed and evolved in the Universe. It has been 
probed observationally since the seminal works of Lilly et al. 
(1996) and Madau et al. (1996, 1998) by using various star- 
formation rate (SFR) indicators, such as luminosities of stel- 
lar continuum at the rest-frame ultraviolet (UV) and nebular 
emission lines (e.g., Ha). Various observational results are 
compiled in Hopkins (2004; H04) and Hopkins & Beacom 
(2006; HB06), in which the cosmology, stellar initial mass 
function (IMF), and dust obscuration correction are unified. 
Their best-fit CSFRD function have been widely utilized not 
only in observational studies (e.g., Karim et al. 201 1) but also 
theoretical studies (e.g., Coward et al. 2008; Tominaga et al. 
2011; Wang etal. 2011). 

The CSFRD at z < 1 has been firmly confirmed by var- 
ious SFR indicators with wide-field survey such as the 
Galaxy Evolution Explorer (GALEX; e.g., Wyder et al.2005; 
Robotham & Driver 201 1). Its characteristic feature is a rapid 
increase with redshift (p* oc (1 + z) 3-4 at z < 1; H04 and ref- 
erences therein). Thus, the CSFRD at z ~ 1 is an order of 
magnitude higher than that in the local universe. By contrast, 
the CSFRD is still uncertain in the higher redshift range (i.e., 
Z > 1), where popular SFR indicators are the rest-frame UV 
continuum stellar emission and IR dust emission. This is sim- 
ply because there are several uncertainties in the estimation of 
the CSFRD from observed data: dust obscuration correction 
for UV continuum, contamination from old stellar population 
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to the IR luminosity, estimation of total IR luminosity, and 
conversion factor from luminosity into SFR. These uncertain- 
ties can result in the well-known inconsistencies between the 
measured statistical quantity and theoretical one by using the 
HB06 CSFRD function such as stellar mass density (SMD, 
p*. e.g., Wilkins et al. 2008; Benson 2010; Choi & Nagamine 
2011) and core-collapse supernova rate (e.g., Horiuchi et al. 
2011). 

The CSFRD has also been calculated theoretically by us- 
ing galaxy formation models: hydrodynamic simulation (e.g., 
Nagamine et al. 2006) or semi-analytic model (e.g., Cole 
et al. 2000; Nagashima & Yoshii 2004; Benson 2010). In 
these models, a galaxy-by-galaxy basis calculation is exe- 
cuted based on a detailed hierarchical structure formation sce- 
nario. Therefore, CSFRD at a certain redshift can be cal- 
culated by simply integrating the SFR of each galaxy at the 
redshift. Although these theoretical models reasonably repro- 
duce the observed luminosity functions (LFs) and the lumi- 
nosity densities (LDs) at the rest-frame wavelength dominated 
by stellar emission, and SMD both at local and high-z, most 
of them underpredict the CSFRD compared to the observa- 
tionally estimated one (e.g., Nagashima & Yoshii 2004; Na- 
gashima et al. 2005; Benson 2010; Lacey et al. 2010). While 
the underprediction of the CSFRD can be attributed to some 
missing unknown physical processes of galaxy formation and 
evolution, it is worth investigating whether or not the uncer- 
tainties in estimating CSFRD from direct observed data can 
be the origin of the disagreement. 

Here, we examine it through a comparison of the observa- 
tional data compiled in H04 with the mock catalog of galax- 
ies generated by one of the semi-analytic models for galaxy 
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formation, so-called the Mitaka model (Nagashima & Yoshii 
2004; see also Nagashima et al. 2005). The Mitaka model 
reproduces a various kind of observations not only for lo- 
cal galaxies including the stellar continuum LFs (Nagashima 
& Yoshii 2004), but also high-z Lyman-break galaxies and 
Lyman-a emitters (Kashikawa et al. 2006; Kobayashi et al. 
2007, 2010). In this work, we focus on the rest-frame UV 
continuum luminosity as a SFR indicator, which is usually 
applied to high-z universe (e.g., Madau et al. 1996, 1998). In- 
vestigation of the other SFR indicators, such as rest-frame IR 
continuum emitted by interstellar dust, will be the subject of 
our future studies. 

As shown in the top and middle panels of Figure Q] the Mi- 
taka model nicely reproduces the measured SMDs and LDs at 
the rest-frame UV wavelengths, which are not corrected with 
interstellar dust extinction, in the redshift range of z = 0-6. 
However, the model prediction to the CSFRD underestimates 
at z > 1 by a factor of ~ 2-3 (i.e., ~ 0.3-0.5 dex) relative to 
the median value of p\ compiled in H04, as shown in the bot- 
tom panel of Fig.Q] It should be emphasized that these H04 
CSFRD data are calculated by using the same observational 
UV LDs plotted in the middle panel of Fig. [TJ which are rea- 
sonably reproduced by our model. Moreover, our CSFRD is 
consistent with the upper limit for given by Strigari et al. 
(2005), which is estimated from an upper limit of the diffuse 
supernova neutrino background with the Super-Kamiokande. 
In this paper, we show that the underestimation of CSFRD can 
be fully attributed to two observational uncertainties, that is, 
dust obscuration correction and SFR conversion used in ob- 
servational studies. Therefore, the underestimation of CSFRD 
does not necessarily imply that theoretical galaxy formation 
models should incorporate some missing unknown physical 
processes. 

This paper is organized as follows. In Section|2] we denote 
the uncertainties in estimating the CSFRD by using observa- 
tional UV LFs. In Section[3] we describe the key prescriptions 
in the Mitaka model which relate to the observational uncer- 
tainties in estimating CSFRD. Then we compare the model 
results with the observed LFs, LDs, and CSFRD in the red- 
shift range of z = - 10 in Section|4] We summarize this work 
and provide some discussions including a new formula for ob- 
scuration correction at the rest-frame UV wavelength in Sec- 
tion |5] Throughout this paper, we adopt the 737 cosmology, 
i.e., H = 70 km s" 1 Mpc" 1 (i.e., h 10 = h/0.1 = 1), fl M = 0.3, 
and Qa = 0.7, and a Salpeter IMF with a mass range of 0.1- 
60 Mq . All magnitudes are expressed in the AB system and 
the wavelengths are given in the rest-frame unless otherwise 
stated. 

2. UNCERTAINTIES IN ESTIMATION OF COSMIC 
STAR FORMATION RATE DENSITY 

In the process of estimating CSFRD at a certain redshift 
from observational data of galaxies, the basic direct observed 
data is their luminosities which can be used as SFR indica- 
tors. First, those are corrected by dust extinction in order to 
obtain intrinsic luminosities, which are expected to more di- 
rectly correlate with SFR. Next the LF is constructed with the 
corrected data which are flux-limited samples and their LD is 
derived by integrating the LF. Lastly, the LD is converted to 
the CSFRD by using the SFR conversion factor. Therefore, 
the uncertainties can be (1) faint-end slope of LF, (2) conver- 
sion factor from luminosity into SFR, and (3) dust obscuration 
correction. 

Here we briefly describe them before we examine the uncer- 



tainties caused by the approaches adopted in H04 and HB06 
to estimate CSFRD from observational data, focusing on the 
UV continuum luminosity as a SFR indicator. 

2.1. Faint-End Slope of Luminosity Function 

In general, the uncertainty of the faint-end slope of LF, a, 
leads that of the integrated LD. This is more significant at 
higher-z simply because only relatively bright galaxies are de- 
tected at higher-z in flux-limited surveys and therefore the un- 
certainty of a becomes larger. 

Figure |2] shows the LDs as a function of luminosity which 
are normalized by the characteristic luminosity L* and num- 
ber density 0* and calculated by using the Schechter function 
form of LF with various a. For a typical value of a ~ -1.5 
in the UV LFs at z = 0-6 (e.g., Oesch et al. 2010), the un- 
certainty of Aa = ±0.3 results in that of the integrated LD to 
log 10 [L/U] --3 with -0.2-0.3 dex (a factor of - 1.6-2). 
This uncertainty can be the origin of the dispersion of LDs at 
each redshift as shown in the middle panel of Fig.Q] 

However, because the UV LD calculated from the Mitaka 
model reasonably lies around the median value of the ob- 
served UV LDs, this is not the origin of the systematic un- 
derestimation of the CSFRD calculated by the Mitaka model 
to that of the H04. Hence, we do not treat this uncertainty 
further in the present paper. 

2.2. Conversion from Stellar Continuum Luminosity into 

SFR 

The intrinsic (i.e., without dust extinction) stellar contin- 
uum luminosity at UV wavelength, L™ l vv , is one of the most 
excellent SFR indicator. This is because it is dominated by 
the stellar radiation from recently-formed massive stars and 
rest-frame UV photons from z = 3 — 10 redshift to optical or 
near-IR in the observer frame at which we can observe most 
efficiently from the ground-based telescopes. 

In the observational studies, the conversion factor from 
L|?' uv to SFR, Csfr, given by Kennicutt (1998; hereafter, 
K98) has been widely utilized: Cf F 9 | = SFR/L^ UV = 1.4 x 
10 _2S M Q yr" 1 (ergs s _1 Hz" 1 )" 1 . H04 and HB06 also adopted 
^sfr- However, it should be noted that, as explicitly stated 
in K98, the linear relation holds only for the galaxies which 
continuously form stars over time scales of 100 Myr or longer 
with the solar metallicity and the Salpeter IMF in the mass 
range of 0.1 - 100 M© (see also Madau et al. 1998). 

Figure[3]shows the time evolution of Csfr at A = 1500 A for 
a constant and exponentially decaying star-formation histories 
(SFHs). It is calculated by using the Schaerer (2003) popula- 
tion synthesis model, which is used in our model to calculate 
UV continuum luminosity of galaxies. As the SFR of a galaxy 
having L^y is evaluated as Csfr x ^vuv trie SFR estimated 
b y c sfr x L l n uv wrongly represents the true SFR. When Cf^ 
is larger than Csfr, the SFR will be overestimated, and vice 
versa. 

In the early phase of a constant star formation, Csfr is sig- 
nificantly larger than Cf^, up to > 1 dex at — 1 Myr age even 
in the case of the solar metallicity Z Q as shown in the top pan- 
els of Fig. [3] This is simply because L™uv continuously grows 
as the number of the massive stars with lifetime of — 100 Myr 
increases before it reaches the equilibrium quantity, which 
corresponds to Cfj^[ for Z = Z . Therefore, adopting Cf^ to 
such young galaxies results in underestimation of SFR for a 
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given ^'"uv The detailed extent of underestimation depends 
on age, metallicity, and wavelength. Conversely, we overes- 
timate the SFR of a galaxy which is old enough to be in the 
equilibrium state and which have sub-solar metallicity up to 
~ 0.2 dex (a factor of ~ 1.6). This can be understood by the 
fact that the color of lower metallicity stars is bluer; smaller 
SFR is enough to produce a certain L'"uv compared to stars 
with larger metallicity. Adopting Cf^ to such evolved galax- 
ies with sub-solar metallicity results in overestimation of SFR, 
and the extent of the overestimation depends on metallicity 
and wavelength. 

In the case of the exponentially decaying SFH, Csfr is not 
constant even at later stage of star formation if the e-folding 
time tsf is short (i.e., tsf < 100 Myr) and progressively de- 
creases with time as shown in the bottom panels of Fig. [3] 
This is because L'"' uv does not decrease as fast with time as 
SFR. Massive stars with lifetime of ~ 100 Myr can contribute 
to L™uv even after the star formation activity quenched in tsf. 
The K98 SFR conversion method for these galaxies where 
SFR decays quickly like starbursts results in either underes- 
timation or overestimation of SFR. 

Csfr should be smaller than Cf^ for higher-z galaxies if 
their UV luminosities reach in equilibrium. This is because 
higher-z galaxies typically have sub-solar metallicity. How- 
ever, if young UV LF is dominated by young galaxies whose 
L™' uv is increasing and/or SFR decays quickly, Csfr should 

be larger than C|^. Therefore, it is not straightforward to 
treat Csfr for high-z universe. We investigate the redshift de- 
pendence of Csfr in this paper in detail by using the Mitaka 
model. 

2.3. Correction of Interstellar Dust Extinction 

In order to utilize the UV luminosity of a galaxy as a SFR 
indicator, it is required to correct its interstellar dust extinc- 
tion. This leads to the most influential uncertainty in the es- 
timation of the CSFRD. This is because the amount of in- 
terstellar dust and its extinction are not easily measurable. 
Moreover, commonly used dust extinction curve rises toward 
the shorter wavelength. The required luminosity correction 
can reach > 0.5 dex at the UV wavelengths. Qualitatively, it 
is expected that higher-z galaxies tend to have smaller dust 
extinction if dust-to-metal mass ratio does not evolve with 
redshift because they have lower metallicity at a fixed stel- 
lar mass (e.g., Mannucci et al. 2010, and references therein). 
However, a realistic extent of dust extinction of high-z galax- 
ies and its redshift evolution are highly uncertain. 

H04 adopted two independent dust obscuration correction 
methods to the various observed LFs: common and SFR- 
dependent obscuration corrections. In both methods, obscura- 
tion correction is executed to observed LFs as statistical quan- 
tity. This means that the observationally bright galaxies would 
be intrinsically brightest. 

In the common obscuration correction, H04 adopted Ay = 
0.52 mag for stars, which is a typical obscuration for the UV 
selected local galaxies (see §2.2 of H04). For the dust ex- 
tinction curve, they adopted the starburst obscuration curve 
given by Calzetti et al. (2000) for all galaxies regardless of 
their luminosities and redshift. Under these assumptions, the 
extinction in magnitude for stellar continuum at the wave- 
length of 1500 A is evaluated as A1500 = 1-33 mag. Appar- 
ently, adopting these quantities can result in an overcorection 
of dust extinction in the case of high-z galaxies which have 



smaller amount of dust than the UV selected galaxies at the 
local universe do. 

In the SFR-dependent correction, obscuration for UV con- 
tinuum at Auv, Auv, of a galaxy with observable UV luminos- 
ity of Lj/.uv is obtained by solving the following transcenden- 
tal equation numerically: 



0.797 logL„, uv + 0.31 8A UV + 0.573 



Auv =X(A uv )log 

Z.OO 

(1) 

where X(1500 A) = 8.935 1 in the case of the Calzetti extinc- 
tion curve. Auv for the SFR-dependent correction as a func- 
tion of observable absolute magnitude at the wavelength of 
Auv, M\jv, is shown in Figure|4]for Auv = 1500 A as well as 
Auv for the common correction. As shown in Figure|4] A 1500 
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■16 mag is adopted to be mag as in H04 although 
the numerical solution for Auv of eq. ([T]) becomes negative. 

This equation is originally derived from an empirical rela- 
tion between £(B- V)g as and Lfir for normal galaxies, blue 
compact galaxies, and starbursts at the local universe, in 
which brighter galaxies have larger obscuration (see refer- 
ences given in §2.2 of H04). However, the distribution of 
these galaxies in the £(6- V)g as -LFiR plane has significant 
scatter around the empirical relation (see Fig. 1 of Hopkins 
et al. 2001). Moreover, the numerical constant in the numer- 
ator in the right-hand side of eq. (UJ) contains the uncertainty 
of SFR/Lfir and SFR/L"' UV . Therefore, this obscuration cor- 
rection methods may lead highly uncertain results at each red- 
shift. Thus, we investigate here the redshift dependence of the 
mean quantity of Auv in detail. 

3. MODEL DESCRIPTION 

In order to examine the cause of the disagreement in the 
CSFRD in theoretical and observational studies, we utilize a 
mock catalog of galaxies generated by a semi-analytic model 
for galaxy formation named "Mitaka model" (Kobayashi et al. 
2007, 2010; updated version of Nagashima & Yoshii 2004). 
The Mitaka model follows the framework of the ACDM 
model of structure formation and calculates the redshift evo- 
lution of the physical quantities of each galaxy at any redshift 
via semi-analytic computation of the merger history of dark 
matter halos and the evolution of baryon components within 
halos. The time evolution of baryon components within halos 
is followed by using physically-motivated phenomenological 
models for radiative cooling, star formation, supernova feed- 
back, chemical enrichment of gas and stars, and galaxy merg- 
ers. As a result, various physical and observational properties 
such as SFR and intrinsic and observable stellar continuum 
luminosities of the galaxies at any given redshift can be ob- 
tained. 

A detailed description of our model is given in Nagashima 
& Yoshii (2004) and Kobayashi et al. (2007, 2010). There- 
fore, here we briefly describe some key prescriptions which 
are closely related to this work. Then, we explain how to 
examine the origin of the difference in CSFRD between ob- 
servation and theory by using the mock catalog of galaxies 
generated by the Mitaka model. 

3.1. Merger Tree of Dark Matter Halo 

1 The definition of X(Auv) here is different from that of Hopkins et al. 
(2001) by a factor of 2.5: X(A UV ) = 2.5Jf Ho P 01 (A uv ). 
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The merging histories of dark matter halos are realized by a 
Monte Carlo method based on the extended Press-Schechter 
formalism (Bond et al. 1991; Bower 1991; Lacey & Cole 
1993). In terms of halo mass function to provide the weight 
for summing merger trees, the Mitaka model adopts the ana- 
lytic functional form given by Yahagi et al. (2004), which is 
a fitting function to their high-resolution iV-body simulation. 

Dark matter halos with circular velocity V c > Vi ow = 
30 km s" 1 are regarded as isolated halos, which corresponds 
to the lower limit of halo mass Mhaio > 10 9 M Q and 10 10 M Q 
at z = 5 and 0, respectively. According to the existence of 
Viow, fainter galaxies (i.e., M1500 > -10 mag) are not well re- 
solved in the Mitaka model. However, this limited resolution 
do not affect the resultant quantities of pisoo and p+ because 
M\soo ~ — 10 mag is much fainter than M* and P1500 well con- 
verges at much brighter magnitudes as shown in the bottom 
panel of Fig. [5] 

3.2. Stellar Continuum Luminosity 

The intrinsic luminosities and colors of galaxies in our 
model are calculated from their SFH and chemical enrichment 
history by using a stellar population synthesis model. We em- 
phasize that we do not assume a simple linear relation among 
the intrinsic luminosity and SFR even at UV wavelengths like 
H04. While the original Mitaka model uses the population 
synthesis model of Kodama & Arimoto (1997), our model 
uses that of Schaerer (2003) to calculate UV stellar contin- 
uum luminosity as in Kobayashi et al. (2007, 2010). Com- 
pared to the Kodama & Arimoto (1997) model, the Schaerer 
(2003) model is a more recent one and covers a wider range 
of metallicity including zero-metallicity. 

It has already found that the revised version of the Mi- 
taka model nicely reproduces all of the observational statis- 
tical quantities of high-z Lyman-a emitters and Lyman Break 
Galaxies (Kobayashi et al. 2007, 2010). 

3.3. Prescription for Dust Extinction 

In order to calculate observational (i.e., with interstellar 
dust extinction) luminosities and colors of model galaxies, 
the optical depth of their internal dust is estimated in our 
model. We make the usual assumption that dust abundance 
of a galaxy is proportional to its gas metallicity and the dust 
optical depth Td is proportional to the column density of met- 
als: 

M C Z C 

T d (X 3—, (2) 

n 

where M c and Z c is the mass and metallicity of cold gas, re- 
spectively, and r e is the effective radius, all of which are ob- 
tained in our model. The normalization of eq. (0 is assumed 
to be a constant and universal regardless of galaxy properties 
and redshift, and has been determined to fit the observed data 
of local galaxies. 

The wavelength dependence of dust extinction r^(A) is 
adopted to be the same as the Galactic extinction curve given 
by Pei (1992). In terms of the dust distribution, our Mitaka 
model assumes the slab dust. Therefore, the amount of ex- 
tinction magnitude A\ for stellar continuum is given by 

10 -o.4A,_ l-exp[-r d (A)] ^ Q) 

Trf(A) 

Note that, while the wavelength dependence of A a calculated 
via eq. (0 is similar to that of Td (i.e., the Galactic extinction 



curve) at Trf <C 1, they are significantly different at tj ^> 1; A\ 
becomes flatter than Td at 7^ ^ 1 as A\ oc log 10 t\. Therefore, 
the wavelength dependence of A\ in our model at Td S> 1 is 
significantly different from that of the Calzetti-law. 

We emphasize that our model calculates star formation and 
chemical enrichment consistently in the framework of hierar- 
chical structure formation. Therefore, both of the time evolu- 
tion of Td of each model galaxy and redshift evolution of its 
mean quantity at a certain redshift can be naturally incorpo- 
rated into our model along with the evolution of M c , Z c and 
r e . 

3.4. Application to the CSFRD Study 

In this work, we calculate the CSFRD in the redshift range 
of z = 0- 10 by using the mock catalog of galaxies generated 
by our model following the prescriptions in H04. We treat the 
mock catalog like the observational data which are compiled 
from literatures and utilized to obtain the CSFRD in H04. We 
emphasize the most striking difference between our catalog 
and the data compiled in H04 is that each model galaxy in our 
catalog has the information of its intrinsic UV luminosity and 
SFR as well as the observable luminosity. 

In terms of obscuration correction, we apply the same 
obscuration corrections with H04 (i.e., common and SFR- 
dependent corrections) into the observable UV LFs at z = 
0-10 given by our model in order to obtain the intrinsic LFs, 
which are represented as common- and SFRdep LFs, respec- 
tively. Their luminosity densities are represented as p^™ and 

P\jy Aep , respectively. We note here that the observable UV 
LF given by our model can be separated into the contribution 
from each model galaxy. Therefore, the intrinsic UV LF can 
be directly provided via simply counting the contribution of 
model galaxies in each magnitude bin. However, we treat the 
observable UV LF as a statistical quantity in order to exam- 
ine how much the H04 obscuration corrections reproduce the 
"true" intrinsic UV LF and p™y. 

The conversion factor of Cfp^ is also compared to that of 
our model in both of each galaxy SFR/L'^uy and all galax- 
ies p+//0™ v . We note again the intrinsic UV luminosity of a 
galaxy is calculated according to its detailed SFH and chem- 
ical enrichment history; therefore, Csfr can vary in every 
galaxy. 

We note here that, while there are several free parameters 
in the phenomenological models for baryon evolution in our 
model, they are already determined to fit the various observa- 
tions of the local galaxies (i.e., B- and /f-band LFs, neutral 
gas mass fraction, and gas mass-to-luminosity ratio as a func- 
tion of B-band luminosity) in Nagashima & Yoshii (2004). As 
their values are kept unchanged in this study, there is no free 
parameter that can be adjusted. 

4. COMPARISON WITH H04 PRESCRIPTIONS 

4. 1 . Dust Obscuration Correction 

We first show the results related to the H04 approaches cor- 
recting interstellar dust obscuration. The intrinsic and observ- 
able UV LFs for A = 1500 A at z = are shown in the top panel 
of Figure [5] The common- and SFRdep LFs and their LDs 
are also shown. Both obscuration correction methods reason- 
ably reproduce the intrinsic UV LF, although there are small 
discrepancies. This might be trivial because both of the H04 
corrections are calibrated by galaxies in the local universe. 
The intrinsic 1500 A LD given by our model is also well re- 
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produced by the H04 methods. 

In order to see the origin of the similarities and discrepan- 
cies between the U V LF and LD given by our model and those 
by the H04 methods, we compare the dust extinction of model 
galaxies as a function of absolute magnitude M1500 with -A 1500 
and Aj™ ep . As shown in Figure [6] the median quantities of 
A1500 for our model galaxies at each magnitude, (A1500), are 
found to take almost constant of (A 1500) ~ 0.6- 1.0 mag in 
-22 mag < M1500 S _ 15 mag. However, model galaxies with 
a certain magnitude have significant dispersion around (A1500) 
up to ~ 1.0 mag. This means that non-negligible number of 
intrinsically bright galaxies lies in observationally faint mag- 
nitudes. The contribution from such observationally-faint and 
intrinsically-bright galaxies to intrinsic LF and LD is com- 
pletely neglected if a single value of (A1500) at a certain mag- 
nitude is adopted as representative extinction for all galaxies 
with the magnitude. Therefore, correcting dust extinction by 
using (Ai5oo)(Mi5oo) results in underestimation of the bright- 
end of intrinsic LF and intrinsic LD. 

Providing the consideration above, we compare the H04 
obscuration corrections with the median quantity (A 1500) of 
the model galaxies. Extinction in the common correction 
of A [goo = 1.33 mag is found to be larger than (A 1500) by 
« 0.3 — 0.7 mag, which are within l-cr error except for the 
brightest magnitude. As shown in the top panel of Fig. [5] 
this over-correction of dust extinction compared to the correc- 
tion with (A i5oo ) is not sufficient to compensate for the con- 
tribution of the observationally-faint and intrinsically-bright 
galaxies at the bright-end of intrinsic LF. However, the over- 
correction at faint magnitudes is enough to cover the under- 
estimation at the bright-end. Hence, the common correction 
results in a similar intrinsic UV LD with that of the Mitaka 
model. On the other hand, extinction in the SFR-dependent 

correction of A [5,^ ) dep is significantly different from (A 1500)- It 
is larger and smaller than (A 1500) in the magnitude range of 
M1500 J; — 18 mag and A/1500 ^ — 18 mag, respectively. The 
over-correction at bright magnitudes leads to the overestima- 
tion of the bright-end of intrinsic LF as shown in Fig. [5] It 
is compensated for the underestimation of number density of 
faint galaxies. As a result, the SFR-dependent correction also 
results in a similar intrinsic UV LD with that of the Mitaka 
model. 

We show the redshift evolution of the intrinsic 1500 A LD 
ratios R\ 500 = Pi5oo/Pi5bo> wnere ' = com or SFRdep, as a 
function of magnitude in Figure [7] As both correction meth- 
ods ignore the contribution of the observationally-faint and 
intrinsically-bright galaxies, the numerical quantities of ^i 50 o 
significantly apart from unity at bright magnitudes. However, 
^j5oo rapidly converge into a certain quantity at faint magni- 
tudes; the asymptotic quantities of R' l500 at z = are found 
to be unity, as seen in Fig. [5] However, it is easily found that 
their asymptotic quantities increase with z and reach ~ 0.3 dex 
and ~ 0.4 dex at z = 3 and z > 4, respectively. In summary, 
the H04 obscuration corrections alone results in the overes- 
timation of CSFRD by ~ 0.1-0.4 dex in the redshift range 
Z = 0— 10 compared to that of our model galaxies. 

4.2. Ratio of UV Luminosity to SFR 
In Figure|8] we show the distribution of model galaxies with 
^1500 — ~15 ma g at z = in Csfr - M500 plane. It is found 
there are two sequences in the plane: a constant sequence 
at ~ CEp^ and widely spread sequence around it with rela- 
tively high and low number densities (> 10~ 4 hj Mpc~ 3 and 



< 10~ 6 hj Mpc" 3 ), respectively. The former and latter se- 
quences correspond to the two distinctive galaxy population 
in our model, that is, quiescently star-forming and starburst 
galaxies, respectively. 

The quiescent galaxies have long star-formation timescales 
(> 1 Gyr) and form disk stars almost constantly in the 
timestep of the Mitaka model. Therefore, their Csfr already 
reach equilibrium quantities, which depend on the stellar 
metallicity. As there is a well-known magnitude-metallicity 
relation (i.e., fainter galaxies have smaller metallicity; e.g., 
Garnett 2002), Csfr decreases toward fainter magnitude. In 
contrast, the starburst galaxies, which starburst activity is trig- 
gered by a major merger of galaxies, have short star-formation 
timescales (tsf < 1 Gyr) determined via dynamical timescale 
of newly formed spheroid. Because their onset times of 
starburst are chosen randomly in the timestep of the Mitaka 
model, they have a variety of Csfr according to the onset time 
and tsf as shown in the bottom panels of Fig. [3] The starburst 
galaxies whose Csfr are larger than Cf^ are relatively young 
and those with Csfr < Cf^R are old enough to have such small 
Csfr- 

We also evaluate the median quantities of Csfr for all of 
model galaxies at each magnitude, (Csfr)- Starburst galaxies 
dominate the bright-end in M1500 at z = (< -20 mag) while 
quiescent galaxies determine (Csfr) except for the bright-end. 
As a combination of the contributions from these two distinc- 
tive populations, (Csfr) is found to be almost constant in all of 
the magnitude range shown in Fig. [8] and slightly smaller than 
Csfr- Therefore, converting L^yy into SFR via C|p^ results 
in overestimation of SFR for most galaxies at z = 0. 

Let us define the effective L™' uv -to-SFR conversion factor, 

C e s i{(z), via: 

qg(z) = pJz)/Pvy(z). (4) 

This is not a conversion factor for each galaxy but the direct 
one from the statistical quantity of intrinsic UV LD p^ v into 
CSFRD p\ by definition. Figure [9] shows the redshift evolu- 
tion of C e s {f R (z) at A = 1500 A. As shown in Figure C^{ 
for A = 1500 A at z = is smaller than Cf^ by w 0.15 dex. 
This result is natural because (Csfr) is smaller than Cf^ at all 
magnitude range of -22 mag < Mi 500 < —15 mag as shown in 
Fig. [8] This implies that converting p^v of the model galaxies 
at z = into p\ by using C|p^ results in overestimation of p\ 
compared to the true quantity by « 0.15 dex. 

The redshift evolution of C e S p R is also shown in Fig. [9] It 
is not simply understandable because of the following effects. 
Toward high redshift, the contribution from young starbursts 
increases because of increasing rate of major merger. This ef- 
fect results in increase of C™ . At the same time, the dynam- 
ical timescale of starburst galaxies becomes shorter at higher 
redshift because of the formed galaxies tend to be smaller in 
size and hence have shorter tsf- This effect results in decrease 
of Cgp{ because Csfr rapidly decaying in the galaxies with 
short tsf- These two competing effects lead to an increase of 

C^FR toward z < 8 and its decrease at z > 8. 

Cgp£ for A = 1500 A is found to be always smaller than 
Csfr by ~ 0.10-0.15 dex. In summary, adopting the usually 
utilized conversion factor of C|p^ for the galaxies at z = 0- 10 
results in overestimation of p\ by s» 0.10-0.15 dex. 
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5. SUMMARY AND DISCUSSIONS 

In this work, we examined the cause of the inconsistency 
in CSFRD between theoretical and observational studies, fo- 
cusing on a SFR indicator for high-z universe, that is, rest- 
frame 1500 A stellar continuum luminosity Lisoq. By using a 
semi-analytic model for galaxy formation so-called the Mi- 
taka model, we found that the underestimation of CSFRD 
seen in theoretical models is originated from the following 
two uncertainties in the process to evaluate CSFRD from ob- 
served 1500 A LD: dust obscuration correction and conver- 
sion from Li5oo to SFR. Uncertainty in the faint-end slope 
of LF is not the origin of the underestimation of CSFRD but 
the origin of dispersion around median CSFRD quantity. The 
methods of obscuration correction adopted in H04 result in 
the overestimation of CSFRD by « 0. 1 -0.4 dex, and the SFR 
conversion used in observational studies also leads to the over- 
estimation of CSFRD by « 0.1-0.2 dex. 

Combined with the fact that theoretical models including 
ours reproduce the observed data of UV LD which is not cor- 
rected by dust extinction and that of SMD, this result indi- 
cates that the inconsistency in CSFRD is not necessarily im- 
plies theoretical models miss some key physical processes in 
galaxy formation. Of course, theoretical models are not per- 
fect yet because there remains observed data to be reproduced 
such as cosmic downsizing (e.g., Cowie et al. 1996). Revis- 
ing theoretical models can be executed through a comparison 
with direct observed data which is not affected by a certain 
model and/or assumption. 

In this section, we provide a brief discussion on the origin of 
the difference in dust extinction at high-z between our model 
and the H04 corrections. Moreover, we present new empirical 
calibrations for dust extinction and SFR conversion as well as 
a recipe how to utilize them in observational studies. 

5.1. Origin of the Difference in Dust Extinction at High 
Redshift 



As described in § 14.11 the H04 obscuration correction meth- 
ods reproduce the intrinsic LD of our model from the observ- 
able UV LF of the model galaxies at z = 0, although they over- 
estimate it at higher redshift. Here we discuss the origin of the 
overestimation at high redshift. 

Our model naturally incorporate the redshift evolution of 
dust abundance because our model calculates chemical en- 
richment of gas in each model galaxy consistently accord- 
ing to its SFH. It can be seen in Figure [TUl which represents 
the redshift evolution of universal effective dust extinction in 
magnitude, A^f^, for our model galaxies defined by: 



A^ Q ( Z ) = -2.5 log [pi 5 oo(z)M n 50o«] 



(5) 



As a consequence of the redshift evolution of metallicity and 
dust abundance, A j{y becomes smaller at higher-z. However, 
such redshift evolution of dust abundance is not incorporated 
into the H04 obscuration corrections. This is the origin of 
their overestimation of the intrinsic LD. 

Similar quantity with has been evaluated from the ob- 
served LD ratio between IR and UV, pir/ puv in the redshift 

range z = 0- 1 (Takeuchi et al. 2005). They found A^q mono- 
tonically increases toward higher redshift, which is the oppo- 
site trend with ours. Examining the origin of the difference is 
interesting but this will be the subject of our future studies. 



5.2. Empirical Calibrations to Obscuration Correction and 
SFR Conversion from UV Luminosity 

Here we propose new empirical formulae which correctly 
correct dust obscuration and to convert from intrinsic UV LD 
Puy into CSFRD pV All of them are derived to reproduce the 
true quantities of the model galaxies and represented by an- 
alytic functions in explicit form. For obscuration correction, 
we show two different formulae; one is for observable LF and 
the other is for its integrated quantity of puv- The former and 
latter formulae are similar to the SFR-dependent and com- 
mon corrections of H04, respectively, but they are described 
to have redshift dependence. 

5.2.1. Conversion from Observable LF into Intrinsic LF 

Let us define Cd us t as an empirical formula to convert ob- 
servable LF into intrinsic one. Cd us t at a certain magnitude 
M1500 is determined to satisfy that the cumulative number den- 
sity of the observable LF of the Mitaka model n obs (< Mi 50 o) 
equals to that of the intrinsic LF « lnt (< M1500- Cdust)- The top 
panel of Figure QT| shows the resultant Cdust as a function of 
M1500 at z = 0. We find that Cdust can be fitted well with the 
following analytic function in the redshift range of z = - 10: 



Cdust(M 150 o; z) = a exp [b\M 1500 -M^, 



15001 



(6) 



where a, b, c, and M" 500 are model parameters and evolve 
with redshift. a and M° 50() have a unit of magnitude, while b 
and c are non-dimensional constant. We adopt a smoothly de- 
clining functional form even for the bright magnitudes where 
there is no model galaxies. Although there are significant dis- 
crepancies between the raw quantities of Cd us t and the evalu- 
ated one via eq. © at M1500 < -23 mag or M1500 > —15 mag, 
the intrinsic UV LF and LD of the model galaxies can be 
nicely reproduced from the observable LF with the analytic 
function of Cdust as shown in Fig.[TT| TableQ]gives the numer- 
ical quantities of the best-fit parameters in the redshift range 
ofz = 0-10. 

The normalization factor of Cdust, a, indicates the maximum 
dust extinction in magnitude at the redshift. a is found to 
gradually increase with redshift toward its peak at z ~ 1-3 
and then decreases as z increases. It is interesting that the 
redshift at which a has its peak roughly corresponds to the 
redshift at which dusty galaxies (e.g., ultra-luminous infrared 
galaxies, sub-mm galaxies, etc.) are mainly found. 

5.2.2. Conversion from Observable UV LD into Intrinsic One 
The conversion factor from observed 1500 A LD pisoo into 



intrinsic one p'^oo * s defined as C { 



dust 



Pi500/P 



1500- 



This re- 



lates to the effective dust extinction A^g defined in eq. Q 



via CtL = dex 



0.4A 



1500 



It is found that C 6 /^ of our model 



galaxies can be well fitted by the following simple analytic 
function with two redshift-independent parameters of a and 
f3 in the redshift range of z = 0- 10 as shown in Fig.fTOl 



q^ t (z) = aexp[-/3(l+z)] + l. 



(7) 



This functional form is motivated by natural expectation that 

Ishift; there is little dust there 
1. The best-fit parameters, 



Cllst approaches unity at high redshift; there is little dust there 



and hence C e / U { = pfWpisoo 
which reproduce p7soo/ /°i500 ratio at z = 0- 10 within ±10%, 
are given in Table [2] 
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While the maximum extinction parameter a in eq. © has its 
peak at z — 1-3, C^ decreases monotonically with increas- 
ing z. The different redshift evolution can be interpreted with 
the fact that the contribution from the galaxies fainter than 
Mj 500 , where the galaxies have a maximum extinction of a at 
the redshift, to intrinsic 1500 A LD is significant; for a typi- 
cal faint-end slope of UV LF a « — 1.4, the contribution from 
the galaxies with L < L„ reaches ~ 0.7 dex as shown in the 
right panel of Fig. [2] Since such faint galaxies have smaller 

Cdust than its peak quantity of a, C e JJ. t progressively decreases 
toward high redshift although a has its peak at z — 1-3. 

5.2.3. Conversion from Intrinsic UV LD into CSFRD 

It is found that the ratio of CSFRD to intrinsic 1500 A 
LD, /di/pj^QQ, can be well fitted with the following simple 
quadratic function in the redshift range of z = 0- 10: 

C&) = Co [ 1 + Q ( 1 + z) + C 2 ( 1 + z) 2 ] , (8) 

where Co, Ci, and C2 are the model parameters and redshift- 
independent constants. Co has the same dimension with 

Cg£{ of M Q yr" 1 (ergs s" 1 Hz" 1 )" , while C x and C 2 are non- 
dimensional constants. With the best-fit quantities given in 
Tab. |2] the analytic function reproduces our model results of 

P*M500 within ±5% - 

We note here that, while the empirical formula given in 
eq. dHJ predicts that C e §J R has a peak at z ~ 6 and progressively 
decreases toward high redshift, it is simply a fitting result and 
does not have any physical moti vatio n. Nevertheless, it might 
be a real trend as discussed in § 14.21 

5.3. Recipe for Converting Observed UV LF into CSFRD 

Here we describe a recipe how to convert observed 1500 A 
LF into CSFRD p\ by using our empirical formulae given in 
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§E2] 

As a first step, intrinsic 1500 A LD p^oo should be calcu- 
lated from the observed UV LF. This can be done by using 
either one of the following two approaches. One approach is 
converting the observed UV LF into intrinsic LF via the em- 
pirical formula of Cdust given in eq. © with the best-fit param- 
eters in Tab. Q] Then, the intrinsic UV LF is integrated over 
magnitude to obtain the intrinsic UV LD, p^oo- The other 
approach is integrating the observed UV LF over magnitude 
first to obtain observable UV LD pisoo, and then it is con- 
verted into (O'^oo by using eq. (|7). Then, one can evaluate p\ 
from pj5 00 with eq. ||8}. 

While the best-fit parameters are provided discretely in red- 
shift for our empirical formulae, it is statistically enough to in- 
terpolate the formulae for a specified redshift. For Cdust given 
in eq. ©, one should execute an interpolation of the param- 
eters to evaluate the adequate Cdust at the desired redshift, in- 
stead of interpolation of Qust itself. 

The numerical quantities of p\, 01500, and Pi^qq of the Mi- 
taka model are compiled in Table Q\ We also present the CS- 
FRD parametric fits to a variety of analytic forms in the liter- 
atures (Cole et al. 2001; Hernquist & Spring el 2003; Yiiksel 
etal. 2008) in Table |U 



We thank Shunsaku Horiuchi and John Beacom for useful 
discussions. The numerical calculation were in part carried 
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tional Astrophysics, CfCA, of National Astronomical Obser- 
vatory of Japan. MARK was supported by the Research Fel- 
lowship for Young Scientists from the Japan Society for the 
Promotion of Science (JSPS). 



REFERENCES 



Benson, A. J. 2012, NewA, 17, 175 

Bond, J. R., Cole, S., Efstathiou, G., & Kaiser, N. 1991, ApJ, 379, 440 

Bower, R. 1 99 1 , MNRAS , 248, 332 

Choi, J.-H., & Nagamine, K. 22012, MNRAS, 419, 1280 

Cole, S. et al. 2001, MNRAS, 326, 255 

Coward, D. M., Guetta, D., Burman, R. R., & Imerito, A. 2008, MNRAS, 
386, 111 

Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 112, 839 
Elbaz, D. et al. 2010, A&A, 518, 29 
Gamett, D. R. 2002, ApJ, 581, 1019 

Hopkins, A. M., Connolly, A. J., Haarsma, D. B., & Cram, L. E. 2001, AJ, 
122, 288 

Hopkins, A. M. 2004, ApJ, 615, 209 (H04) 

Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 (HB06) 

Horiuchi, S., Beacom, J. F., Kochanek, C. S., Prieto, J. L., Stanek, K. Z., & 

Thompson T. A. 201 1, ApJ, 738, 154 
Karim, A. et al. 201 1, ApJ, 730, 61 
Kennicutt, R. C. Jr. 1998, ARA&A, 36, 189 (K98) 

Kistler, M. D., Yiiksel, H., Beacom, J. E, Hopkins, A. M., & Wyithe, J. S. B. 

2009, ApJ, 705, 104 
Kennicutt, R. C. Jr 1998, ARA&A, 36, 189 

Kobayashi, A. R. M., Totani, T., & Nagashima, M. 2007, ApJ, 670, 919 
Kobayashi, A. R. M., Totani, T, & Nagashima, M. 2010, ApJ, 708, 1119 
Kodama, T. & Arimoto, N. 1997, A&A, 320, 41 
Lacey, C. G., & Cole, S. 1993, MNRAS, 262, 627 

Lacey, C. G., Baugh, C. M., Frenk, C. M., & Benson, A. J. 201 1, MNRAS, 
412, 1828 

Lilly, S. J., Le Fevre, O., Hammer, F, & Crampton, D. 1996, ApJ, 460, LI 



Madau, R, Ferguson, H. C, Dickinson, M. E., Giavalisco, M., Steidel, C. C, 

& Fruchter, A. 1996, MNRAS, 283, 1388 
Madau, P., Pozzetti, L., & Dickinson, M. 1998, ApJ, 498, 106 
Mannucci, F, Cresci, G., Maiolino, R., Marconi, A., & Gnerucci, A. 2010, 

MNRAS, 408, 2115 
Murphy, E. J., Chary, R.-R., Alexander, D. M., Dickinson, M., Magnelli, B., 

Morrison, G., Pope, A., & Teplitz, H. I. 2009, ApJ, 698, 1380 

Murphy, E. J., Chary, R.-R., Dickinson, M., Pope, A., Frayer, D. T, & Lin, 

L. 2011, ApJ, 732, 126 
Nagashima, M., & Yoshii, Y. 2004, ApJ, 610, 23 

Nagashima, M., Yahagi, H., Enoki, M., Yoshii, Y, & Gouda, N. 2005, ApJ, 
634, 26 

Oesch, P. A., et al. 2010, ApJ, 725, L150 
Pei, Y. C. 1992, ApJ, 395, 130 

Robotham, A. S. G., & Driver, S. P. 2011, MNRAS, 413, 2570 
Strigari, L. E., Beacom, J. F, Walker, T. P., & Zhang, P. 2005, JCAP, 4, 17 
Takeuchi, T. T., Buat, V., & Burgarella, D. 2005, A&A, 440, L17 
Tominaga, N., Morokuma, T., Blinnikov, S. I., Baklanov, P., Sorokina, E. I., 

& Nomoto, K. 2011, ApJS, 193, 20 
Wang, F. Y, & Dai, Z. G. 201 1, ApJ, 727, 34 

Wilkins, S. M., Trentham, N, & Hopkins, A. M. 2008, MNRAS, 385, 687 
Wyder, T. K., et al. 2005, ApJ, 619, L15 
Yahagi, H., Nagashima, M., & Yoshii, Y. 2004, ApJ, 605, 709 
Yiiksel, H., Kistler, M. D., Beacom, J. E, & Hopkins, A. M. 2008, ApJ, 683, 
L5 



8 



Kobayashi, Inoue, & Inoue 



TABLE 1 

Fitting Parameters for Our Formula to 
Correct Dust Extinction, Cd ust 



redshift z 


a [mag] 


b 


c 


M 1500 [ ma §] 





1.847 


-0.1139 


1.566 


-19.49 


1 


2.466 


-0.1181 


1.281 


-22.65 


2 


2.344 


-0.2625 


1.086 


-21.86 


3 


2.235 


-0.2844 


1.184 


-21.91 


4 


2.032 


-0.2115 


1.428 


-22.15 


5 


1.822 


-0.1971 


1.599 


-22.16 


6 


1.583 


-0.1673 


1.737 


-22.14 


7 


1.496 


-0.1630 


1.800 


-22.08 


8 


1.415 


-0.1900 


1.641 


-22.16 


9 


1.228 


-0.1239 


1.895 


-22.39 


10 


1.109 


-0.1898 


1.634 


-22.17 



TABLE 2 

Fitting Parameters in Our Formula of CfIJ a and (fJJL 



r*SS 

'-dust 




r*ff 

SFR 


C o [M /yr/(ergs/s/Hz)] 


a 13 




c, 


2.983 0.3056 


-5.915 x 10" 5 


7.294 x 10- 4 


-28.01 



TABLE 3 

LD AND CSFRD OF THE MlTAKA MODEL 



redshift z 




P1500 


^1500 


P'*/pl500 




[h 70 M yr" 1 MpiT 3 ] 


[h 70 ergs s~' Hz -1 Mpc -3 ] 


[h 70 ergs s 1 Hz 1 Mpc 3 ] 


[M© yr 1 (ergs s" 1 Hz-'r 1 ] 





2.599 x 10- 2 


8.090 x 10 25 


2.616 X 10 26 


3.213 x 10" 28 


1 


6.825 x 10- 2 


2.531 x 10 26 


6.593 x 10 26 


2.697 x Kr 28 


2 


8.285 x 10- 2 


3.665 x 10 26 


7.910 x 10 26 


2.261 x Kr 28 


3 


6.828 x 10~ 2 


3.498 x 10 26 


6.396 x 10 26 


1.952 x 10" 28 


4 


4.926 x 10~ 2 


2.752 x 10 26 


4.518 x 10 26 


1.790 x 10" 28 


5 


3.124 x 10- 2 


1.902 x 10 26 


2.818 x 10 26 


1.642 x 10" 28 


6 


1.759 x 10~ 2 


1.119 x 10 26 


1.540 x 10 26 


1.572 x Kr 28 


7 


9.105 x 10~ 3 


5.988 x 10 25 


7.849 x 10 25 


1.521 x 10 -28 


8 


4.524 x 10~ 3 


3.129 x 10 25 


3.865 x 10 25 


1.446 x 10 -28 


9 


2.069 x 10~ 3 


1.641 x 10 25 


1.898 x 10 25 


1.261 x 10" 28 


10 


9.514 x 10- 4 


8.245 x 10 24 


9.150 x 10 24 


1.154 x 10- 28 



TABLE 4 

CSFRD PARAMETRIC FITS TO THE FORM OF THE FOLLOWING LITERATURES 



Ref. Functional Form parameter value 



Cole etal. (2001) p\(z) = (a+bz)h[\ + (z/c) d ] a 0.0389 

b 0.0545 

c 2.973 

d 3.655 

Hernquist & Springel (2003) p» (zj = p X 2 / 1 1 + a( X - 1 ) 3 exp (j3 X ' /4 )] /i, 0.030 

X(z)=[H(z)/H ] 2 P a 0.323 

13 0.051 

Yuksel et al. (2008) p\{z) = po[(l +z) arl + {(1 +z)/B}^ + {(1 +z)/C}-f>]Vv t p 0.0258 

B = (l+ Zi y- a ^, C = (1+z 1 )C 3 - q Vt(1+z 2 ) 1 ^t a 1.6 

P -1.2 

7 -5.7 

Zl 1.7 

Z2 5.0 

r] -1.62 
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1 2 3 4 5 6 7 



red shift z 

FIG. 1 . — Redshift evolution of SMD (top), rest-frame UV (i.e., A = 1500-2800 A) LD (middle), and CSFRD (bottom). The solid curves in each panel are the 
predictions of the Mitaka model. The thin solid curve in the top panel is the SMD for the model galaxies with M„ > 10 8 Mq. The dashed curve in the bottom 
panel are the best-fitting Cole et al. (2001) functional form to the H04 p\. The symbols with error bars are observational data. The observational data shown in 
the top panel is compiled in Yabe et al. (2009) and that in the middle and bottom panels is done in H04, respectively. In the bottom panel, the boxes (circles) are 
evaluated by using a common (SFR-dependent) obscuration correction in H04 from the same observed data plotted in the middle panel. 
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-3 -2.5 -2 -1.5 -1 -0.5 -3 -2.5 -2 -1.5 -1 -0.5 

log 10 [L/U] log 10 [L/U] 

FIG. 2. — LD contribution from the galaxies with brighter than the horizontal axis of L/L* calculated by using the Schechter function with various faint-end 
slope of LF a. As represented by the legends, the solid, long-dashed, short-dashed, dotted, and dash-dotted curves are a = —1.2, -1.4, -1.6, -1.8, and -1.9, 
respectively. Top and bottom panels are the same, but with a vertical axis of linear and log scales, respectively. The left and right panels are the normalized LD 
by 0„L* and the total LD, respectively. 
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SFRoc e xp[-t/T SF ],Z = Z 
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FIG. 3. — Time evolution of Csfr = SFR/L'"™ calculated by using the population synthesis model of Schaerer (2003) with a Salpeter IMF in the mass 
range of 0. 1 — 60 Mq . Top: Csfr for constant star formation. The solid, long-dashed, and short-dotted curves are Csfr in the case of the stellar metallicities of 
Z/Zq = 1 , 1 /5 , and 1 / 50, respectively. Bottom: Csfr for exponentially decaying star formation with e-folding time of tsf in the case of Z = Zq . The solid, long- 
dashed, short-dashed, long dash dotted, and short dash dotted curves are Csfr at log 10 [rsF/yr] = 10, 9, 8, 7, and 6, respectively. The horizontal solid line in all 



panels represents the most popular conversion factor utilized in observational studies given by Kennicutt (1998), C™ = 



1.4 X 10~ 28 M Q yr 1 (ergs s~' Hz~ 
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-26 -24 -22 -20 -18 -16 

M 150 o-5log 10 h 70 [ABmag] 



FIG. 4. — Interstellar dust extinction in magnitude at A = 1500 A, A1500, as a function of observable (i.e., with dust extinction) absolute magnitude M1500 
adopted in H04. The thick curve is A^qq in the case of the SFR-dependent obscuration corrections, which is obtained by solving the transcendental equation of 
eq. (TJ numerically. The thin horizontal line represents A 1500 m me case of the common obscuration correction. A 1500 at M1500 > — 16 mag for the SFR-dependent 
correction is mag by definition as in H04. 




FIG. 5. — 1500 A LFs (top) and LDs (bottom) as a function of absolute magnitude at z = 0. The thin solid- and doited curves are the intrinsic and observable of 
the Mitaka model, respectively. The thick solid- and dashed curves are the common- and SFRdep-LFs and LDs, which are calculated with applying the common 
and SFR-dependent obscuration correction of H04 to the observable LF of the Mitaka model. 
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log 10 [n / (h^ Mpc" 3 )] 




-22 -20 -18 -16 

M 1500 - 5log 10 h 70 [ABmag] 

FIG. 6. — Distribution of the model galaxies with < -15 mag at z = in the A1500-M1500 plane. The color scale shows the color-coded values of the 

number density of the model galaxies in unit area. The open-circles with error bars connected with lines are the mean and lcr for the model galaxies in each 
magnitude bin. The obscuration corrections of H04 are also shown in solid curve and dotted line for the common and SFR-dependent, respectively. The dashed 
lines indicate the contours of MJj' 0() corresponding to -22, -20, and -18 ABmag from left to right. 
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FIG. 7.— The 
Z = — 7 from top 



LD ratios of Rjj™p SFRdcp = Pi5oo SFRdep / Pi'soo as a function of absolute magnitude for the wavelength of A = 1500 A in the redshift range of 
-left to bottom-right. The line styles are the same as in Fig. [5] 
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FIG. 8. — Same as Figure[6]but for distribution in the Csfr-Misoo plane. C|™ = 1.4 X 10~ 28 Mq yr~' (ergs s _I Hz~') 1 is also shown as dotted horizontal 
lines. 
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FIG. 9. — Ratio of CSFRD and intrinsic 1500 A LD of the Mitaka model as a function of redshift. The open circles are the model results and the solid curve 
shows the best-fit quadratic function for them defined by eq. (8), which best-fit parameters are compiled in Tab.|2] Note that the dynamic range of the vertical 
axis (0.3 dex) is much smaller than that of Fig.[8](2.0 dex). 
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FIG. 10. — Same as Figure[9]but for the ratio of the observable 1500 A LD to intrinsic one of the Mitaka model as a function of redshift. The parameters for 
the best-fit analytic function defined by eq. (7) and represented as thick solid curve are also written in Tab. [2] 
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FIG. 1 1 . — Top, middle, and bottom panels show Cdust, LF, and ratio of the LD calculated from the Mitaka observable LF -fCdust to the Mitaka intrinsic LD, 
respectively, at z = for the wavelength of A = 1500 A. In the top panel, the filled circles are the numerical data calculated from the intrinsic and observable 
1500 A LFs of the Mitaka model, while the dashed curve is analytical fit by using eq. |6} with the numerical quantities in Tab.[T] Cdust = at Muoo < -22 mag is 
reflected by the fact that there is no model galaxy with such bright observable magnitudes. In the middle panel, the solid- and dotted curves are the intrinsic and 
observable LFs of the Mitaka model, respectively, while the dashed curve is the Mitaka observable LF +Cd us t- 



